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Abstract 

We study the antiferromagnetic RP 2 model in four dimensions. 
We find a second order transition with two order parameters, one fer- 
romagnetic and the other antiferromagnetic. The antiferromagnetic 
sector has mean-field critical exponents and a renormalized coupling 
which goes to zero in the continuum limit. The exponents of the fer- 
romagnetic channel are not the mean-field ones, but the difference can 
be interpreted as logarithmic corrections. We perform a detailed anal- 
ysis of these corrections and conclude the triviality of the continuum 
limit of this model. 
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1 Introduction 



The non-perturbative formulation of non-asymptotically free, interacting field 
theories in four dimensions is yet to be accomplished. The conventional anal- 
ysis, for A</> 4 and O(N) theories, yields triviality in four dimensions [p]]. That 
is, once the continuum limit is taken, correlation functions factorice as the 
Wick's theorem prescribes for the gaussian theory. A possible way, in or- 
der to obtain a model with a non-trivial continuum limit, is to introduce 
antiferromagnetism (AFM). Gallavotti and Rivasseau have considered AFM 
actions to change the ultraviolet limit of 4 theories 0. From a Statistical 
Physics point of view, a great variety of AFM models in three dimensions has 
been studied to obtain different qualitative behaviour from that of the cor- 
responding ferromagnetic (FM) models ||. In four dimensions, recent works 
have studied the possibility of new universality classes if AFM is added [|J . 

The AFM RP 2 model has recently been studied in three dimensions, 
due to its exotic properties [|, For instance, it has a disordered, unfrus- 
trated ground state. Even more, it seems to present a full breaking of the 
action's 0(3) symmetry ||. Perturbative studies of this Spontaneous Sym- 
metry Breaking (SSB) pattern, yield the 0(4) universality class |§. If this 
prediction holds true in four dimensions, the fate of the model is triviality. 

However, this theoretical prediction has been questioned in three dimen- 
sions by Monte Carlo (MC) simulations ||. Therefore, the study of the 
triviality of this model in four dimensions is very interesting. It also would 
help to enlighten the situation in three dimensions. We will see, although, 
that a detailed analysis of the MC simulation of this model indicates the 
triviality of its continuum limit through the appearance of logarithmic cor- 
rections to the divergences of the observables of the theory. A special form of 
the Finite Size Scaling (FSS) analysis which includes logarithmic corrections 
will be used to deal with these corrections. 

We define the model and observables in section §|2|, where we also describe 
the techniques we have used to measure the critical exponents. The results 
of the MC simulation are presented in section §|3]. The model exhibits a 
phase transition at a negative coupling with two independent order param- 
eters. One of this channels, the staggered one, presents mean-field critical 
exponents, but the other, ferromagnetic, presents deviations. We show in 
sections §[3] and §^ how the discrepancies with the mean-field behaviour can 
be interpreted as logarithmic corrections. 
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2 The model 



We shall consider the RP 2 = S 2 /Z2 (real projective space) spin model in 
four dimensions. Our basic variable is a three component normalized spin 
Vi, interacting through a gauge Z2-invariant action. As a local symmetry 
cannot be broken, these are effectively RP 2 variables (only the direction of 
the vectors is relevant). We consider a hypercubic lattice, with first neighbour 
coupling: 



<ij> 




The ferromagnetic (positive coupling) model presents a first order transi- 
tion at (3 ~ 0.94. The ground state consists of spins parallel or antiparallel 
to an arbitrary direction, and the SSB is SO(3)/SO(2). The analysis of the 
antiferromagnetic counterpart is trickier, given the more complicated nature 
of the ground state. Let us call a lattice site, labeled by (x,y,z,t), even 
or odd according to the parity of x + y + z + t. In the ground state, ev- 
ery even/odd spin is parallel or antiparallel to an arbitrary direction, while 
odd/even spins lie randomly on the perpendicular plane. The corresponding 
SSB is SO(3)/SO(2), which calls for the 0(3) universality class. However, 
fluctuations induce an interaction between spins on the randomly-plane sub- 
lattice. In three dimensions, this seems enough to break the remaining 0(2) 
symmetry ||, and to change the Universality Class. 

2.1 Definition of observables 

The natural RP 2 variables are given by the traceless tensorial field Tj 

Tf = vfv? - ±fi°f>, (2) 
whose lattice Fourier transform will be represented by T: 

T p = J2 e M~ip-r) T r . (3) 

r 

We will work in a L 4 lattice with periodic boundary conditions. We define 
two order parameters, according to the discussion of the ground state above, 
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the intensive staggered (ferromagnetic) magnetization, as the sums of tensors 
on even sites minus (plus) those on odd sites, or equivalently 



1 - 1 - 

M s = -^T (7ri7ri7rj7r ) , (M = yjT( 0) o,o,o))- (4) 

As no spontaneous symmetry breaking can occur on a finite lattice, in a MC 
simulation one needs to measure 0(3)-invariant operators. For the magneti- 
zation and the susceptibility, we define 

M = (VtrM 2 ) , x = £ 4 (trM 2 ), (5) 



and analogously with the staggered observables. 

A very useful quantity for a triviality study is the Binder cumulant. For 
this model, we define 



2 V5 (trM 2 )' 

which, in the infinite volume limit, becomes 1 in the broken phase and in 
the symmetric one. The cumulant for the staggered magnetization is defined 
analogously. 

Another very interesting quantity is the second momentum correlation 
length defined as U 

U = (w(tt/L)) ' (7) 

where F is the mean value of the trace of T squared at minimal momentum 
(27r/L in any of the four directions). For £ s we use Xs and F B , analogously 
defined from T at momentum (2ir/L + 7r, tt, tt, it) and permutations. 

The field theoretical definition of the renormalized coupling constant can 
now be introduced 

g R = V M (L/Z L ) d , (8) 

where d is the dimension of the lattice. We will consider the renormalized 
couplings associated with the two different sectors. 

In addition, we measure the energy, which is needed for the spectral 



density method flQfl, invaluable for extrapolating MC measures to a neigh- 



bourhood of the critical coupling. 
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2.2 Standard Finite Size Scaling 

To study critical exponents, we have used a method specially suited to the 
measurements of anomalous dimensions ||, [ff|. Let us consider the mean 
value of an operator O, measured in a size L lattice, at a coupling value (3 
in the critical region. Let t be the reduced temperature {(3 — /3 C )//3 C . The 
standard FSS formula states that |TT] 

{0{L,t)) = {0{t))To{s{L,t)) , 8 (L,t) = jjpr, (9) 

where 0(t) means 0(oo,t) and To is a smooth function. We suppose that 
the values of L and are large, so that we ignore scaling corrections in 
Now, we have (0(t)) ~ t~ x ° , which is the definition of the critical exponent 
Xo, and ~ t~ u , and we can write L x °l v = (0(t))s x °/ u . This allows to 
write @ as 

(0(L,t)) = L x °l»g {s). (10) 
Applying (|10"D to the correlation length, it gives £,(L,t) = LQ^(s), so that 

x(L,t) = Qj£ = gt( 8 ), (n) 

and s = (x). From eq. (|T0|) we have 

(0(L,t)) = L*ol»g (g^(x)) = L x °/»f (x), (12) 
and we finish up with the useful expression 

(0(L,t))=r/7 ^j+.., (13) 

where the dots stand for possible scaling-corrections. Let us denote 

(Q(rL,t)) 

Qo ~ (0(L,t))' (M) 

we can produce with it a sensible measure of critical exponents: 

Qo\ Q , =r = r*°/» + ■ ■ ■ . (15) 

Therefore, from simulations of lattice sizes L and rL, we can extract the 
critical exponent xq/v from the quotient (|H]), measured at the point where 
one correlation length is r times the other. 
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L 



r(Xs) 



Measures (xlO 3 r (Xs)) 



4 
6 

8 

10 
12 
16 
20 
24 



0.644(16) 
1.189(20) 

2.26(3) 
3.65(11) 
5.53(17) 

10.7(2) 



30 

60 

64 

19 

21 

39 

5/5 

4/2 



17.0(10)/2.31(11) 
25. 5(10)/3. 53(17) 



Table 1: Total number of measures and the corresponding integrated cor- 
relation times r for Xs- in the larger lattices the data of the overrelaxed 
simulations (right) are separated from those of Metropolis (left) by slashes. 



3 The Monte Carlo simulation 

To simulate the system, we have used a standard three-hits Metropolis algo- 
rithm, with an uncorrelated change proposal, achieving approximately a 50% 
of acceptance. The lattice sizes have been L = 4,6, 8, 10, 12, 16, 20 and 24. 
For the larger sizes, 20 and 24, we have combined Metropolis with an over- 
relaxed update, described in Appendix §[A], to decrease the autocorrelation 
time. The overrelaxed algorithm is not able to decrease the dynamic critical 
exponent z, but nevertheless we save total CPU time when compared with 
the simple Metropolis simulation. 

The runs have been distributed over several Workstations. We display 
in Table [l] the integrated autocorrelation time for Xs and the number of 
measurements performed for every lattice size. Every two measurements 
are separated by 10 sweeps, each consisting of either one Metropolis update 
or one Metropolis plus three overrelaxed updates when we use the latter 
algorithm. 

3.1 About the order parameters 

The RP 2 model presents a second order phase transition at (3 ~ —1.34. The 
ferromagnetic and staggered magnetizations defined in equation (]5|) are zero 



6 



below the transition. To show that they are real order parameters, we should 
ensure that they do not vanish in the broken phase when V — > oo. In Fig. p] 
we plot the values of M s and M at /3 — —1.5 for the lattice sizes L = 8, 12 and 
16. It is clear that both magnetizations reach an asymptotic value different 
from zero in the thermodynamical limit in the broken phase. 

Z 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




0.005 0.010 0.015 

1/L 2 

Figure 1: Asymptotic values of M (straight line) and M s /10 (dashed line) 
from their values for L — 8, 12 and 16 at /3 — —1.5. 



3.2 Critical exponents 

We have calculated the critical exponents v and rj for the two different chan- 
nels using (fI5|), which yields x x = 7 for the susceptibility, and xm = ~P for 
the magnetization. To calculate u, we use xq^/q^ — v + 1. All along this 
paper we shall take r = 2. 

We obtain the anomalous dimension rj through the scaling relations: 

(2 - V )u = 7 , 2/3 = u(d-2 + V ). (16) 
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The resulting values for the 77 exponent from these two relations, will be 
denoted by rj x and r\ M respectively. 

As far as the exponent v is concerned, we expect the same critical expo- 
nent for both correlation lengths, the ferromagnetic £ FM and the staggered 
£. We have found that the measures for £ are more accurate so we have used 
this variable as correlation length. 

We plot in Fig. |2| an example of how this method works in both channels. 
Notice that Q M i takes the value 2"^ s ^ u ~ d when Q £ = 2. 

(s) 
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Figure 2: Quotients of M s 2 and M 2 as a function of the quotient of £(L). The 
horizontal straight lines correspond to mean-field behaviour. The symbol 
sizes are proportional to the lattice sizes. 



The resulting exponents are shown in Table [2]. After the name of the 
exponents their mean- field values are shown in square brackets 0. The high 
accuracy reached on the measures of the 77 exponents, is due in part to the 
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Stag 


gered 


Ferromaj 


metic 


Li, Lt2 


i/[0.5] 


V X [0] 


Vm[0] 


V X [2] 


Vm[2] 


4,8 


0.527(8) 


0.009(4) 


0.008(4) 


1.959(8) 


1.963(9) 


6,12 


0.524(6) 


0.006(3) 


0.008(3) 


1.976(5) 


1.981(6) 


8,16 


0.512(4) 


0.008(2) 


0.005(2) 


1.968(3) 


1.973(3) 


10,20 


0.491(7) 


0.013(4) 


0.001(4) 


1.964(6) 


1.969(6) 


12,24 


0.496(9) 


0.007(5) 


0.009(5) 


1.960(8) 


1.965(7) 



Table 2: Estimations for the critical exponents of the AFM RP 2 model. 



strong statistical correlation between and Qm 2 - 

We obtain a value for v compatible with the mean-field prediction as 
well as for the magnetic exponents of the staggered sector. However, in the 
ferromagnetic channel, our exponents are close but not compatible with those 
given by mean-field theory. 

Another possible interpretation of these values is that they could be effec- 
tive critical exponents due to the presence of strong logarithmic corrections 
in the FM sector. To check this, let us suppose that there are logarithmic 
corrections only in the susceptibility. We do not take into account here the 
fact of possible logarithmic corrections to £ as we use the values of the quo- 
tients measured at r = 2 value. This is an approximation that holds for large 
L. We address to section §^ for a more complete treatment of the logarithmic 
corrections. So that 

X ~L m (lnLr, (17) 

Qx= X^ = 2 m fML\ m ^ 2 m h ™ (lg) 



Xl \ hi L 

where m = 7/V. The effective exponents obtained with the standard FSS 
m' L = lnQ x /ln2, can be written as 

m T = m + m . 19 

L In 2 K J 

The fit discarding the 4,8 pair for the ferromagnetic susceptibility gives 

m = 0.09(9), m = -0.13(7), x 2 / d -°- f - = 0.04/2. (20) 
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The fit (PO]) indicates that the exponents of the ferromagnetic channel in Ta- 
ble |2| are compatible with a mean-field situation with logarithmic corrections 
in the susceptibility. 



3.3 Critical temperature 

The thermodynamical critical temperature of our system can be estimated 
from the crossing points of the Binder cumulants for the different lattice 
sizes [12|. To obtain /3 c (oo), we can extrapolate according to the formula 

&(oo) - p c (L) « L~ l /\ (21) 

It should be noted that (|2T| ) is only a first approximation, because it 
does not include any logarithmic corrections. In section § [4.3| we will be able 
to measure /3 c (oo) taking into account logarithmic corrections. However, it 
hardly modifies the value of the critical coupling obtained with this method. 

We show in Fig. ^ the Binder cumulants for the staggered and ferromag- 
netic channels. 

We have fitted the crossing points of the Binder cumulants of the L = 8 
lattice with lattices L = 12, 16, 20 and 24. The results of the fit are 

Staggered Ferromagnetic 

/3 c (oo) = -1.3426(3) /3 c (oo) = -1.3421(6) 
X 2 /d.o.f. = 1.2/2 xVd.o.f. = 1.0/2 

Both values are compatible, and, as we expect one transition point, we 
take the value of /3 c (oo) with lower error, that is 

/3 c (oo) = -1.3426(3). (22) 

Considering the crossing of the Binder cumulants of the L — 10 lattice 
with the larger lattices scarcely change the numbers. 



4 Logarithmic corrections of the RP 2 model 

As we have previously shown, the values obtained for the critical exponents 
for the four dimensional AFM RP 2 model, are compatible with those pre- 
dicted by mean-field plus logarithmic corrections. They appear more clearly 
in the FM channel because in this sector there is no power-law divergence of 
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the susceptibility, while in the staggered sector the logarithmic divergence is 
added to a power-law one. We need a modification of the standard FSS to 
include these corrections. 

4.1 FSS with logarithmic corrections 

Let us consider an observable 0(t) whose behaviour near the critical point 
would have a logarithmic contribution 

(0{t)) ~r x \lntf. (23) 
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We will follow [^] to take into account the logarithmic corrections: the 
scaling variable s(L,t) of @ is now substituted by £(L, 0)/£(t): 

(0(L,t)) = (0(t))J- (^^). (24) 

Formula (|24| ) coincides with (^) below the upper critical dimension, where 
£(L, 0) ~ L, otherwise it can also take into account the logarithmic correc- 
tions to the finite volume correlation length. 
Let us suppose that, to leading order, 

£(L,0)~L A (lnLf, (25) 

where a and (3 are two exponents that depend on the theory. For the O(N) 



models, a — 1, (3 — j [11]. The transition at finite L takes place when t is 



such that f (L, 0) ~ £(t). If f (t) ~ t~ u \ lnt| p , then, employing (g5|), 

t~ L- &/u (\nL) {D -P )/v . (26) 
Below four dimensions a = 1 and (3 = v = 0. Now, we have 

>(L>t)s f^.%^, (27) 

so that making use of fl2"6|), and with a change of variable similar to ([TT]) , we 
obtain 

(0(L,t)> = L->(lnL)^^-^ + Vo (^4^) > ( 28 ) 

\L a (lnL)^/ 

which is the equation analogue to (|IB|). We follow now the same method as 
in the standard FSS case: we compute the quotient 

Q = = ^ r i + ( ^ )+ * . (29) 



(0(L,t)) V InL/ fnM^ 

Measuring Qo at the point ti where 

t(L,t L ) ~ 2 ^ 



(30) 
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with h L = 1 + we find 

Q Q (t L ) = 2 & */»h x L /u0 - 9)+s . (31) 

This is the new expression that substitutes ( 13 ) when there are logarith- 
mic corrections, from it we can extract the exponent xjv. We no longer have 
to measure Qo where the quotient of correlation lengths is 2, but instead 
where it equals 2 a bP L . 



4.2 Staggered channel: renormalized four-point 
coupling 

We proceed now to calculate the renormalized four-point coupling of our 
theory. The limit we are interested in is 

g R = Van g R (L, f3 c (oo)). (32) 

The evolution of g R at the critical temperature with L is shown in Table 
[| where we have used the value for /3 c (oo) of (p2j). 



L Staggered Ferromagnetic 

~8 3.16 (3) -0.19 (3) 

10 2.84 (4) -0.19 (5) 

12 2.61 (5) -0.21 (5) 

16 2.34 (4) -0.23 (5) 

20 2.08 (9) -0.23 (11) 

24 1.95 (13) -0.15 (17) 



Table 3: The renormalized four-point coupling g R (L, /3 c (oo)). 



When hyperscaling is violated by logarithms, g R 
apply the FSS formula (121 



\nt\-P HI. Let us 



9n(L,t) = gR^J 7 , 



HP, 



L & (\nLf 



(33) 
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The scaling behaviour with L at the critical point is 



We can eliminate the t dependence in ([34j) if lim^ Fg^z) ~ z~^' v , which 
makes the | lnt|~^ factor disappear. After that, we employ relation (|26| ) and 
get finally 

g K {L,0)^(hxL)-P, (35) 

which is also directly obtained from ( pq) . We can obtain p fitting the values 
of Table |3| to the functional form (|35|). The fit for all lattice sizes, L > 8 
yields (Fig. 1) 

p = 1.07(6), x7 d -0- f - = 0.8/4. (36) 

The result ( |36"D implies triviality for the staggered sector of the AFM RP 2 
model in four dimensions. The renormalized coupling goes to zero because 
of logarithmic corrections, exactly in the same way as in the ferromagnetic 
0(N) models, for which p — 1 [T5[ . 



The behaviour of the ferromagnetic channel is however rather different 
from that of the staggered sector. From the data of Table we cannot 
conclude an asymptotic value for the renormalized coupling. We have seen 
above that the logarithmic corrections are very strong in this channel, and 
maybe for that reason, the renormalized coupling is very hard to measure. 
To conclude about the triviality of this sector, we will try to study in full 
detail the logarithmic corrections of this channel, following the FSS analysis 



derived in section 84.1 



4.3 Computation of the logarithmic corrections 
4.3.1 Correlation lengths 

To parameterize the logarithmic corrections, we need a lot of, in principle, 
unknown exponents: a, /3, v, . . . Therefore it will be necessary to make a few 
assumptions about these exponents. 

The most important exponents are a and (3, (|25|) , because we have to 
measure at the points which satisfy (|30|). Assuming a mean-field plus log- 
arithmic corrections scenario, we expect a — 1. To find out the value of 
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exponent j3, we use the prediction (|25| ) for the correlation length in a finite 
lattice at /3 c (oo): 

ln^ = C + /31n(lnL). (37) 

We have performed the fit for both correlation lengths, the ferromagnetic 
(£ FM ), and the staggered (£) one. In order to monitorize subleading effects, 
we have compared the fits with L > 8, and L > 10. We have found that 
the fit for £ FM is more stable with growing lattice sizes. In Table f|, we 
show the fit parameters. The infinite volume critical coupling, and the fit- 
parameter errors have been estimated from the increment in one unit of the x 2 
function. Comparing with the previous determination of the critical coupling 
(22), both determinations are consistent and of similar accuracy, although 



logarithmic corrections to scaling were not considered previously. Our value 
for the exponent (3 is consistent with the predicted value for the O(N) models, 
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Fit 


X 2 /d.oi. 


P 


&(oo) 




-^min = 8 


2.1/3 


0.21(2) 


-1.3423(3) 


-^min = 10 


0.3/2 


0.16(3) 


-1.3424(6) 






0.8/3 


0.22(4) 


-1.3425(3) 


-^min = 10 


0.2/2 


0.17(8) 


-1.3424(3) 



Table 4: Fits for the logarithmic corrections to the correlation lengths at the 
critical point. 



(3 = 0.25, specially the obtained from the ferromagnetic correlation length. 
4.3.2 Magnetic operators 

We shall try to control other logarithmic corrections by making use of the 



results of section § [4.1| . From (|29| ) we know that 

Q x = 2 & ^hl /u($ - p)+ \ (38) 



when measuring the quotient at the point which verifies the condition fl30|). 
The logarithmic corrections are given in terms of the unknown exponents 
v and 7, but, as we assume mean-field exponents, 7 = 0, v — 0.5, we can 
reduce ( j38|) to 

lnQ x = 7 Info. (39) 
In a similar way, if we take the magnetization, 

Qm = 2- & V»hZ< 3/v0 -° )+0 , (40) 

or (mean-field: (3=1) 

\ n Q M = -2 In 2 + Rlnh L , (41) 

where R = —2(/3 — v) + (3. Therefore, from every pair of lattices L, 2L, we 
can obtain the exponents 7, R. This is shown in Table [S]. 

Notice that the 7 and R values are very close to their corresponding values 
for the magnetization in 0(3) and 0(4) models. It should be understood that 
these values (7 ~ 0.5 and R ~ 0.25) are calculated for the order parameter 
on the fundamental (vectorial) representation of the O(N) group. However 
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0.6 0.8 

ln(ln L) 



1.2 



Figure 5: Determination of the exponent f3 of the FSS formulas from the 
behaviour of £ FM and £ measured at the mean values of f3 c from the fits with 
L > 8. 



the critical exponents of the FM magnetization in the RP 2 model and those 
of the order parameter in the tensorial representation for the 0(N) models 
are the same at the mean-field level. Let us remark that the 7 ~ 2R just 
means that (M 2 ) ~ (M) 2 . 

This final result from the FM sector completes the conclusion that we 
obtained after examining the renormalized coupling of the staggered sector 
■O: the RP 2 model is trivial due to the logarithmic corrections to the 



m 



mean-field behaviour. 

The question of the SSB pattern remains unsolved as the ferromagnetic 
susceptibility is only logarithmically divergent. We recall that the power-law 
behaviour was crucial to check the symmetry breaking in three dimensions pf . 
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7 



K 



4,8 
6,12 
8,16 
10,20 
12,24 



0.45 (2) 
0.45 (1) 
0.49 (1) 
0.53 (2) 
0.52 (3) 



0.23 (6) 
0.22 (6) 
0.24 (3) 
0.26 (1) 
0.26 (2) 



Table 5: Exponents of the logarithmic corrections of the FM sector of the 
RP 2 model. 

5 Conclusions 

We have examined the triviality question of the four dimensional AFM RP 2 
model, which presents a second order transition. A very interesting feature 
of this model is that it presents two different order parameters. A detailed 
study of these two sectors reveals that the model has a trivial continuum 
limit. We have been able to calculate explicitly the logarithmic corrections 
to the mean-field behaviour by means of a FSS analysis also valid at the 
critical dimension of the model. 
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A The overrelaxed algorithm 

Overrelaxation is a local microcanonical update algorithm. It makes the 
maximum change in the variable at a given point without modifying the 
energy. We will now describe how this method works in our model. It is easy 
to see that 




(42) 
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so that the change T, — > T[ must be such that the trTjN is conserved, 
where N is the sum of the tensors at the neighbour points of i. If we set 
T' = RTR -1 , the conservation of the energy is then expressed by 

[R, N] = 0. (43) 

Besides, we must ensure that the new tensor belongs to RP 2 , so that the 
change in T is associated with a change in v: v' = Rv. As the vectors are 
normalized, this puts also the condition of unitarity on the matrix R . 

In order to fulfill these two conditions, let us write N as N = UAU~ l , 
where A is the matrix of eigenvalues of N, and define C = U~ l RU . Then, 
the updating conditions, written in terms of the matrix C are 

[C,A] = , C- 1 = C + . (44) 



As A is a diagonal matrix (|44|) implies that C has to be also diagonal, and 
C 2 = 1, which means that its three eigenvalues will be ±1. We have reduced 
our updating process to a choice of the matrix C. Here enters the second 
characteristic of the overrelaxed algorithm: the change in the vector v should 
be maximum, which can be achieved by minimizing the value of the squared 
scalar product 

A=(v- v'f = (v ■ Rv) 2 = (v ■ Cvf, (45) 

where 

v = U^v = (xi, x 2 , x 3 ). (46) 

To do the update, we then have to take the three numbers c, = ±1 that 
minimize the quantity 

A = (cix 2 + c 2 x\ + C3S3) 2 . (47) 

To sum up, the overrelaxed algorithm consists of calculating the matrix N 
of nearest neighbours, its eigenvectors to obtain U, and then looking for the 
minimum of the combinations in (fT7). 



It is easy to see that this algorithm verifies detailed balance: if we make 
v — ► v' — > v", so that v" = UC'v' and v' = Cv, we have to minimize 

(«' ■ C'v') 2 = (v ■ CC'Cv) 2 , (48) 

and the last expression was minimized by C, so that C = CC'C or C = C. 
Therefore v" = v. 
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